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We show that the correlations of electrons with a fixed energy in metals have very anomalous 
time and space dependences. Due to soft modes that exist in any Fermi liquid, combined with the 
incomplete screening of the Coulomb interaction at finite frequencies, the correlations in 2-d systems 
grow with time as t 2 . In the presence of disorder, the spatial correlations grow as the distance 
squared. Similar, but in general weaker, effects are present in 3-d systems and in the absence of 
quenched disorder. We propose ways to experimentally measure these anomalous correlations. 
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Equilibrium time-correlation functions are an essen¬ 
tial concept in statistical mechanics [T] . They describe 
the spontaneous fluctuations of a system in equilibrium, 
and together with the partition function they provide 
a complete description of the equilibrium state. Via the 
fluctuation-dissipation theorem they also describe the lin¬ 
ear response of the system to external fields, and they are 
directly measurable by means of scattering experiments. 

An old, and seemingly plausible, assumption is that 
microscopic correlations decay on time scales that are 
much faster than macroscopic observation times. Var¬ 
ious concepts depend on this assumption, for instance, 
the notion that the BBGKY hierarchy of classical ki¬ 
netic equations can be truncated 2]. An analogous as¬ 
sumption underlies the Kadanoff-Baym scheme of deriv¬ 
ing and solving quantum kinetic equations and its gener¬ 
alizations [31 HI- The assumption of a separation of time 
scales is also important in other areas, e.g., in signal pro¬ 
cessing 03 M- For time-correlation functions it implies 
that they decay exponentially for large times. Equiva¬ 
lently, their Laplace transform is an analytic function of 
the complex frequency z at z = 0. The discovery of the 
non-exponential decay known as long-time tails (LTTs) 
urn, and the related breakdown of a virial expansion 
for transport coefficients mm thus came as a consid¬ 
erable surprise [T2], since it showed that the assumption 
is in general not true. Rather, many time-correlation 
functions decay only algebraically, i.e., they have no in¬ 
trinsic time scale. This scale invariance is reminiscent of 
the behavior of correlation functions at critical points; 
however, it occurs in entire phases, as opposed to just 
at isolated points in the phase diagram, and therefore 
is referred to as ‘generic scale invariance ’ [HHE]. The 
underlying physical reason is either conservation laws, or 
Goldstone modes that lead to a slow decay of some long- 
wavelength fluctuations and, via mode-mode-coupling ef¬ 
fects, affect the decay of other degrees of freedom. An 
example is the shear stress in a classical fluid, which is 
not conserved, yet its time-correlation function decays al¬ 
gebraically as l/t d / 2 for long times t in a d-dimensional 


fluid since it couples to the transverse momentum, which 
is conserved. As a result, the Green-Kubo expressions 
for various transport coefficients diverge in dimensions 
d < 2, and the hydrodynamic equations become nonlocal 
in time and space; for a review, see Ref. llbl 

In classical systems in equilibrium, LTT effects, while 
qualitatively very important, are rather small quantita¬ 
tively and become pronounced only at times so large that 
the correlation function is already very small overall. In 
non-equilibrium classical systems the effects are much 
more important m m- In equilibrium quantum sys¬ 
tems the corresponding effects can also be much larger, 
especially in systems with quenched disorder, where the 
quantum LTTs are often referred to as “weak-localization 
effects” [I1]II8]. Still, the correlation functions consid¬ 
ered to date decay as functions of time, albeit more slowly 
than a separation-of-time-scales argument would suggest. 
In this Letter we show that in a quantum system as sim¬ 
ple as interacting electrons with no quenched disorder, 
i.e., the simplest model of a metal, there are correlations 
that not only do not decay exponentially, but actually 
grow with time, and in some cases also with distance. 
This surprising result is a consequence of generic soft, 
or slowly decaying, excitations in a Fermi liquid in con¬ 
junction with the incomplete screening of the Coulomb 
interaction at nonzero frequencies. It is a dramatic il¬ 
lustration of the fact that the impossibility of separating 
microscopic and macroscopic time scales, which is present 
in classical kinetics, holds a fortiori in quantum systems. 

In quantum statistical mechanics it is useful to con¬ 
sider correlation functions that depend on one or more 
imaginary-time variables r € [0,1/T], with T being the 
temperature, or on the corresponding imaginary Mat- 
subara frequencies, icj n = 2i.TrT(n + 1/2) for fermions, 
and id n = 2iTrTn for bosons (n integer). Functions de¬ 
fined for imaginary Matsubara frequencies can be ana¬ 
lytically continued to all complex frequencies, and the 
underlying real-time dependence can be obtained by an 
inverse Laplace transform. The observables in a fermion 
systems can be expressed in terms of expectation val- 
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ues of products of field operators A(x,t) and xf)(x,r) 
that depend on the position x in addition to r. Spin 
is not essential for our purposes, and we suppress it 
for now. Let us consider binary products of rjv and "0, 
and an imaginary-time Wigner operator W ( X , x ; T, r) = 
+ x/2,T +t/2)^{X -x,T- t/2). In a field- 
theoretic formulation, ^(x,t) and ^(a:,r) correspond 
one-to-one to fermionic (i.e., Grassmann-valued) fields 
^}(x,t) and ip( x , T ) [15], in terms of which we define a 
Wigner field 

W(X,x;T,t) = ${X + x/2 1 T+t/2) 

xip(X — s/2, T — t/2) (1) 

in analogy to the operator W. In common applications 
of real-time Wigner operators or fields, X and T cor¬ 
respond to the “average” or “macroscopic” (presumed 
to be slow) length and time scale, and x and r to the 
“relative” or “microscopic” (assumed to be fast) scales. 
The definition of the Wigner field reflects the assumption 
that it is possible and useful to separate these two scales 
B GO- In terms of it, the fluctuating particle number 
density is given by n(X,T) = W(X,T',x = 0,r = 0), 
and the equilibrium single-particle Green function by 
G(x,t) = (W(X,T',x,t)), where (...) denotes an av¬ 
erage taken with the action governing the fermion sys¬ 
tem. If the average is taken in a non-equilibrium state, 
(W) also depends on X and T. In a real-time formal¬ 
ism, with macroscopic time T and microscopic time t, 
( W(X,x;T,t )) = —iG < (X,x;T,t) is the Green func¬ 
tion G < defined in Ref. [3]. Its Fourier transform with 
respect to the microscopic variables, g < (X,T-p,ui), is 
often interpreted as the density of particles with mo¬ 
mentum p and energy ui at the space-time point (X,T) 
HHS]. Switching back to imaginary time and frequency, 
this identifies 

r i/T 

p{X, iuj n ) =T dT dr e^" 1 " W (X, T,x = 0, r) 

Jo 

= Y,M X )J>nA X ) (2) 

a 

as the density of particles with energy uj n at point X, 
i.e., a spatial energy distribution, p is related to the 
full phase-space distribution by summing over all mo¬ 
menta and averaging over the macroscopic (imaginary) 
time. Restoring spin, its spatial Fourier transform reads 
in terms of fermion fields 

p(k,iu n ) = j dX e~ ik X p{X,iu n ) 

= ^2' t P n ^(P + k /2)fp n ,a{p- k/2) , (3) 

P,(T 

where 4>n(p) = \fT]V f dx e~ tpx i/>(x, t) and ip n (P ) = 
y/rJV f dx e lpx il>(x, r), with px = p ■ x — w„t, J dx = 


f dx f 0 1/T dr , and V the system volume, p depends on a 
macroscopic wave vector k, but a microscopic frequency 
uj n . This is in contrast to the number density n, which de¬ 
pends on two macroscopic variables. To verify the physi¬ 
cal interpretation of p we note that its expectation value 
determines the density of states N(w) via 

N(u>) — — — Im(p(fc = 0, iuj n — > w + i0)) (4a) 

7 r V 

The zeroth frequency moment gives the particle number 
N. With r/ = 0 + the usual convergence factor |2T] and 
np(uj) the fermion distribution function we have 


N = Te iu)nT] {p{k = 0, ioj n )) = V j dun?(u)N(u) , 

n 

(4b) 

and the first frequency moment gives the energy E car¬ 
ried by the particles [22] . 

E = TY / e i “ nV iUn(p(k = 0,iu Jn )) 

n 

= V J duj tif(w) u N (w) . (4c) 


Let us now consider the four-fermion correlation func¬ 
tion C PP (X - Y ; iu> n ,icvm) = (8p{X,iu} n )8p{Y,iu m )), 
with Sp = p — (p). This is motivated by two consider¬ 
ations. First, in light of the above interpretation of p, 
C P p provides information about the correlations of en¬ 
ergy levels in the Fermi system: It is the second moment 
of the energy density distribution. Second, the quantity 
i/(k,iu! n ) = p(k,iu! n ) — p(k,—iuj n ) can be interpreted, 
in a technically precise sense, as an order parameter 
(OP) for the Fermi liquid [2H1 [24]. C vv {k\iix n ,iuj rn ) = 
(8u(k, iuj n ) 8v{—k, iu m )) is thus the (longitudinal) OP 
susceptibility in an ordered phase. We will come back 
to this interpretation below. Writing C pp in imaginary- 
time space, and using time translational invariance, one 
easily sees that it consists of two distinct contributions. 
One piece (which one would call “disconnected” in a di¬ 
agrammatic representation) is proportional to 8 nm , and 
a second, “connected” one, is proportional to T [25]. We 
focus on the connected piece by putting ca m = — Lo n and 
eliminate the trivial factor of temperature by defining 

C(k, iix n ) = ^ (8p(k, iuj n ) 8p(—k, -iuj n )) , (5) 

which has a well-defined zero-tenrperature limit [26]. 
Since C depends on a microscopic time scale, the 
separation-of-time-scales assumption would suggest that 
the analytic continuation C(k,iuj n —> z ) is an analytic 
function of the complex frequency z at z = 0, correspond¬ 
ing to exponential decay in imaginary or real frequency 
space. From ordinary LTT physics one might expect that 
a coupling between the fast and slow degrees of freedom 
will lead instead to a nonanalytic function of the form z a , 
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which would lead to a LTT of the form l/t a+1 . We find 
that neither of these expectations is correct in general: 
In a Fermi liquid with a Coulomb interaction in d = 2 
the real-time dependence of the correlation function C is 


p, n-m 



p-k y n-m 


C(k —► 0, t) oc k 2 ln(n/\k\) t 2 (d = 2) (6a) 


where n is the screening wave number. That is, C in¬ 
creases with time as t 2 ; i.e., the correlations get stronger 
with increasing time. This behavior is cut off by a 
nonzero wave number k or, equivalently, by a finite sys¬ 
tem size L; it is valid for times t <C \J Lk/vfk (and t 
much larger than the microscopic time scale), with vf 
the Fermi velocity. In d = 3 the behavior is a LTT, 

C(k —► 0, t) oc 1/np t (d = 3), (6b) 


which is valid for t <C L/v f- For asymptotically large 
times C decays exponentially, with the rate of decay go¬ 
ing to zero as the wave number approaches zero or the 
system size goes to infinity. 

Also of interest are the spatial correlations. For ioj n —> 
0, i.e., for particles close to the Fermi surface, the spatial 
correlations decay only algebraically, 


C(x, iuj n — > 0) oc 


y/k/vp y/r 
k/v | r 2 


(d = 2) 
(d = 3) 


(7) 


for distances r = \x\ 1/k. 

These results hold for clean systems. In the presence 
of quenched disorder, the effects are even stronger. The 
time dependence in d = 2 is the same as in the clean 
case and given by Eq. (6a), but in d = 3 the correlation 
function does not decay with time for t <C L 2 /D, 


C(k — ► 0,i) oc 1/D 2 |fc| (d = 3), (8) 


where D is the diffusion coefficient that characterizes the 
diffusive electron dynamics. The spatial correlations for 
particles near the Fermi surface grow quadratically with 
the distance and remain constant in d = 2 and d = 3, 
respectively, 


C(x, ioj n —► 0) oc 


(1/D 3 ) r 2 

1/D 3 


(d = 2) 
(d= 3) . 


(9) 


These expressions are again valid for distances large com¬ 
pared to the microscopic length. Note that for particles 
at the Fermi surface (ita n = 0) the spatial correlations at 
large distances in d = 3 get cut off only by the system 
size. In d = 2 they grow as the square of the distance for 
distances less than the localization length or the system 
size, whichever is smaller. 

We now sketch the origin and derivation of these sur¬ 
prising results, and then discuss their physical signifi¬ 
cance as well as ways to check them experimentally. We 


FIG. 1: Diagrammatic representation of the correlation func¬ 
tion C(k,iu n ) for clean systems within (a) the effective field 
theory of Ref. |23j . and (b) many-body perturbation theory. 
In (b), solid and double-dashed lines denote electronic Green 
functions and dynamically screened Coulomb potentials, re¬ 
spectively. Notice that no frequency is transferred at the ex¬ 
ternal vertices (heavy dots); this reflects the fact that this is 
not a contribution to the usual density correlation function. 
The frequency conservation at the internal vertices is as usual. 


first consider clean systems. The correlation function C, 
Eq. (§, can be calculated in various ways. In the frame¬ 
work of the effective held theory developed in Ref. [23] the 
leading contribution is given by the one-loop diagram 
shown in Fig. 0 a). The advantage of this framework is 
that the renormalization-group analysis of the effective 
theory guarantees that the result is the leading behav¬ 
ior; higher-loop diagrams will change the prefactor, but 
not the functional form of the result. Alternatively, the 
same result can be obtained from many-body perturba¬ 
tion theory m via the diagram shown in Fig. 0 b ); how¬ 
ever, there is no such guarantee within that formalism. A 
simplified analytical expression for either diagram, which 
has the correct scaling behavior, at T —> 0 is 


C(k, ioj n ) oc / dpp' 


d—l 


dto 


(to 2 + Vpp 2 ) 2 


(U( P ,to)) 2 . 

( 10 ) 


Here A is an ultraviolet momentum cutoff, and U(p , to) 
is the dynamically screened Coulomb interaction. The 
factor of 1 /(u 2 + v^p 2 ) 2 represents the soft fermionic 
modes. The strongly singular behavior discussed above 
results from a combination of the latter and the incom¬ 
plete screening of the Coulomb interactions at nonzero 
frequencies. A short-range interaction still leads to sin¬ 
gularities, but they are weaker than in the Coulomb 
case; the corresponding behavior is obtained by replac¬ 
ing U(p,u>) in Eq. (10) by a constant. The limit on the 
time regime where Eq. (6a) is valid results from the most 
singular behavior of C in d = 2 being restricted to fre¬ 
quencies oj n larger than the plasma frequency. Note that 
the latter can be made arbitrarily small by going to small 
wave numbers (or to large system sizes at k = 0). 

For disordered systems, an appropriate effective field 
theory is the generalized nonlinear sigma model that has 
been studied extensively in the context of metal-insulator 
transitions [281430] . The relevant one-loop diagram is still 
given by Fig. 0a), but the nature of the propagators 
is diffusive rather than ballistic. Within the framework 
of many-body perturbation theory the diagram shown 
in Fig. 0b) needs to be dressed with diffusion poles in 
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elaborate ways. With any calculation method the net 
result is that the factor l/(w 2 + v^p 2 ) 2 in Eq. (10) gets 
replaced by a diffusion pole to the fourth power, and the 
dynamically screened Coulomb potential gets modified to 
reflect the diffusive nature of the electron dynamics. 

In summary, we have shown that a correlation func¬ 
tion that describes particle-number fluctuations with a 
fixed energy in a Fermi liquid is a very singular function 
of space and time. In the zero-wave-number limit, and 
in a 2 -d system, the correlations grow quadratically with 
time for times small compared to the inverse plasma fre¬ 
quency, and they decay only as 1/t in d = 3. For particles 
near the Fermi surface, the correlations decay only alge¬ 
braically with distance in clean systems, Eq. ([T]). In the 
presence of quenched disorder, they grow quadratically 
with distance in d = 2 and remain constant in d = 3, 
Eq. §). Remarkably, this behavior is cut off only by 
the system size. In what follows we add some discussion 
remarks to put this remarkable behavior in context. 


(1) As mentioned after Eq. (4c), the correlation func¬ 


tion C can be interpreted as an OP susceptibility for the 
Fermi liquid. An interesting analogy in this context is the 
corresponding OP susceptibility in a classical Heisenberg 
ferromagnet. Due to a coupling between the longitudi¬ 
nal and transverse magnetization fluctuations the longi¬ 
tudinal magnetic susceptibility xl (be., the OP suscep¬ 
tibility) for 2 < d < 4 diverges everywhere in the or¬ 
dered phase as l/k 4 ~ d 13T1 . This results from a one-loop 
contribution to Xl that is a wave-number convolution of 
two Goldstone modes, each of which scales as an inverse 
wave number squared. Diagrammatically this contribu¬ 
tion has the same form as Fig. |lja) . To see the origin 
of the stronger effects discussed here, consider the spa¬ 
tial variation of xl as a function of the distance. Setting 
all wave-number components except for k x equal to zero, 
k = (k x , 0,...), we have 


Xl(M 


°° ) = /' 


dk x e l * x xiXk x ) oc \x 


3-d 


(ii) 


That is, the correlations grow with distance for 2 < d < 
3. Our results for the Fermi-liquid OP susceptibility C 
are in direct analogy to this result if one makes the fol¬ 
lowing adjustments: (i) For the time or frequency depen¬ 
dence, replace k x by the frequency and put k = 0. (ii) 
Realize that the relevant propagator in the quantum field 
theory (35] scales as a soft mode squared, see Eq. (10). In 


the many-body calculation, this is apparent from the tri¬ 
angular fermion loops in Fig. 0b) , each of which scales as 
a ballistic propagator squared, (iii) Take into account the 
incomplete screening of the Coulomb interaction, which 
enhances the effect compared to the naive expectation 
that the quantum result should correspond to the classi¬ 
cal one in an effective dimension d e s = d + 1. 

(2) The temporal and spatial dependences of the cor¬ 
relation function C are quite different: The underlying 
correlation is a function of two points in space, but four 


points in time; translational invariance implies that one 
and three of these, respectively, are independent. The 
time dependence of the function C results from having 
integrated over two of the three independent time argu¬ 
ments, which is justified by the physical interpretation 
of the function C. A related point is that we study the 
behavior of C for both frequency arguments approach¬ 
ing the Fermi surface, w n = — oj m —1 0, rather than for 
large frequency differences. The spatial dependence of 
C, on the other hand, has the same structure as in usual 
two-point correlation functions. An important result is 
that the spatial correlations become more and more long 
ranged as the Fermi surface is approached. 

(3) In the classical-magnet analog the strong fluctu¬ 
ations eventually lead to an instability of the ordered 
phase at the ferromagnetic transition. In the present 
case, this suggests the possibility of a transition from a 
Fermi liquid to a non-Fermi liquid with a vanishing den¬ 
sity of states at the Fermi surface [32]. However, whether 
such a transition is actually realized in any given system 
is a question the current theory cannot answer. 

(4) We suggest two ways to experimentally observe the 
effects discussed here: (i) A direct measurement of the 
energy density distribution, e.g., in a cold-atom system, 
and (ii) measurements of the distribution of the local 
density of states, e.g., by tunneling, which is related to 
the correlation function C as can be seen from Eq. (4a I. 

(5) Studies of the distribution of the local density of 
states in disordered metals [33] (33] have calculated a dif¬ 
ferent correlation function:, viz., the disorder average of 
the disconnected piece of the correlation function C UL , de¬ 
fined after Eq. (4c). The effects considered were thus en¬ 
tirely determined by disorder fluctuations and vanish in 
the clean limit. In contrast, the effects considered here 
are caused by the electron-electron interaction, and some 
of them are further enhanced by disorder. 
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